Theory of biphoton generation in a single-resonant optical parametric oscillator 

far below threshold 
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We present a quantum-theoretical treatment of biphoton generation in single-resonant type-II 
parametric down- conversion. The nonlinear medium is continuously pumped and is placed inside a 
cavity which is resonant for the signal field, but nonresonant for the idler deflected by an intra-cavity 
polarizing beam splitter. The intensity of the classical pump is assumed to be sufficiently low in 
order to yield a biphoton production rate that is small compared to the cavity loss rate. Explicit 
expressions are derived for the rate of biphoton generation and for the biphoton wave function. The 
output spectra of the signal and idler field are determined, as well as the second-order signal-idler 
cross-correlation function which is shown to be asymmetric with respect to the time delay. Due to 
frequency entanglement in the signal-idler photon pair, the idler spectrum is found to reveal the 
longitudinal mode structure of the cavity, even though the idler field is not resonant. 

PACS numbers: 03.67.-a,42.50.-p,42.50.Ar,42.65.Lm 



I. INTRODUCTION 

In parametric down-conversion, a pump photon of fre- 
quency ujp incident on a medium with a second-order 
nonlinear susceptibility \ is split into two photons with 
lower frequency [lj. Spontaneous parametric down con- 
version produces photon pairs that can be entangled in 
many degrees of freedom [2J. The resulting two-photon 
state, consisting of a signal photon at frequency u) s and 
an idler photon at frequeny uji — uj p — lu s , is often called 
a biphoton. Upon post-selection on the idler photons, 
the process provides a source of heralded single photons 
which represent the principal resource in many quantum 
information processing protocols like quantum cryptog- 
raphy [H, 0, [BJ or linear optics quantum computation Q . 
Quantum networks have been proposed 0, [H that rely on 
stationary atoms or ions as information processing nodes 
and on single photons to transmit information via optical 
fibers. First building blocks of this scheme have already 
been realized [9|, [ltj . For an efficient atom-photon cou- 
pling, the photon bandwidth has to match the linewidth 
of the atomic transition which is by orders of magnitude 
smaller than the bandwidth of the photons emitted in 
spontaneous parametric down-conversion. 

In order to reduce the photon bandwidth, cavity- 
enhanced parametric down-conversion can be applied. A 
bright source of heralded narrow-band single photons was 
experimentally realized [ll[ using a double-resonant op- 
tical parametric oscillator (OPO). The nonlinear crystal 
was placed inside a cavity resonant for both the signal 
and the idler field, and in the output a single narrow- 
band longitudinal signal mode was selected with an ex- 
ternal frequency filter. Clearly, to conditionally achieve 
single-photon generation, the OPO has to be operated in 
the regime far below threshold where the production rate 
of down-converted photons is small compared to the loss 
rate of the cavity. A number of preceding experiments for 
biphoton generation in a double-resonant OPO far below 
threshold have been performed [H [U, [H, Gj| EE El > 



and the first theoretical description was given in Ref. [13| . 
Recently, the theory of the double-resonant OPO has 
been extended by analyzing the conditionally prepared 
single-photon state [HI], and by providing a multi-mode 
treatment which is valid for both pulsed and stationary 
pump fields 

To ensure reliable operation of a quantum network, 
continuous photon emission over a long period of time is 
essential. For this purpose, an active stabilization of the 
OPO is necessary which proves to be a complicated task 
in the double-resonant case while it is easier to achieve for 
a single-resonant cavity. Continuous biphoton generation 
in a single-resonant OPO far below threshold has recently 
been demonstrated in our group [20( using a setup where 
only the signal mode experiences resonance enhancement 
in the cavity while the orthogonally polarized idler mode 
is nonresonant due to deflection by an intra-cavity polar- 
izing beam splitter. Additional passive filtering with the 
help of an external cavity can select a single longitudinal 
mode and will thus enable the generation of narrow-band 
single photons. 

To our knowledge, a theoretical treatment of the single- 
resonant OPO far below threshold has not been per- 
formed so far. The present paper aims to fill this gap. 
Based on the concepts of the pioneering theoretical stud- 
ies of spontaneous parametric down-conversion 

nana, 

we provide the theoretical background for our experi- 
mental results [2(J. Different from the approaches used 
for the theoretical description of the double-resonant 
OPO [HElEi, free - field quantization of the idler field 
is inevitable for our scheme. 

This paper is organized as follows: In Sec. II, the 
basic equations for describing the nonlinear interaction 
between the quantized signal and idler fields are pro- 
vided. The biphoton production rate and the biphoton 
wave functions are derived in Sec. Ill by applying the 
standard perturbative treatment in the Schrodinger pic- 
ture Q. The results are used in Sec. IV to derive the 
spectral properties of the emitted radiation and to study 
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the the second-order signal-idler cross-correlation func- 
tion. Sec. V concludes the paper by establishing the con- 
nection to real experimental situations. 



II. BASIC EQUATIONS 

A. The interaction Hamiltonian 

We consider type-II parametric down-conversion in 
a nonlinear crystal of length I that is pumped by a 
monochromatic linearly polarized classical field of fre- 
quency up. The crystal is assumed to be placed in- 
side a cavity which is resonant for the signal field, but 
not for the idler field, polarized orthogonal to the signal. 
Fig. 1 shows a schematic picture of the corresponding 
experimental setup [20]. The central frequencies u>s and 



nonlinear 
crystal 




FIG. 1: (Color online) Scheme of the considered experiment. 
The cavity acts as a one-sided resonator for the signal field, 
while the idler field is non-resonant due to deflection at an 
intra-cavity polarizing beam splitter. The signal-idler cross- 
correlation function is determined by delayed coincidence de- 
tection. 

u>i of the signal and idler field depend on the properties 
of the birefringent nonlinear crystal and are determined 
by energy and momentum conservation known as phase- 
matching [l[ 



UJp = lj s + wi, 
k p (ujp) = k a (u>s) + ki(uJi). 



(2.1) 
(2.2) 



p , fv s , and ki are the wave vectors of the pump, 



Here, kp, k 

signal, and idler waves which are collincar in the consid- 
ered setup, with the signal leaving the cavity in positive 
x-direction. The standing-wave fields of the pump and 
the signal inside the cavity are composed of two compo- 
nents, propagating in negative and positive x-direction, 
respectively. Because of the phase-matching conditions, 
only those components contribute to the parametric in- 
teraction that have the same propagation direction as the 
idler wave. Neglecting the vector notation of the fields, 
the positive- frequency part of the electric field in the rel- 
evant component of the classical pump inside the crystal 
can be written as 



Ep,cr(x,t) = E P e 



i[k p (cjp) x—Ljpt] 



(2.3) 



The weak signal and idler fields inside the crystal are 
described by the operators E 



(+) 

S,cr 



4^ and Ey; 



( + ) 



Ej , respectively, denoting the positive- frequency part 
of their co-propagating field components. The interaction 
Hamiltonian can be written in the simplified form 



r i _ X 

Hint. — 777 



(+) 



— / dx(Ep, cr E Scr E Icr + E Pcr E Scr E Icr ), 

where the second-order nonlinear susceptibility \ is 
frequency-dependent [l|. Before utilizing Eq. (|2.4p . we 

need to find explicit expressions for the operators E s cr 

and Ej 



B. The free-field operators 

For later use, we start by providing the operators for 
the signal and idler fields in free space. The positive- 
frequency electric field operator of a wave with transverse 
cross-section A propagating freely in x-direction is given 

by E(+\x,t) 



limi^oo £\ o 



rGUe 



"*) with 



j=U y 2£ LA" 1 ^ 

ujj = 2njc/ L where L is the quantization length and <Zj 
denotes the photon annihilation operator of mode j. Let 
us consider the signal field and introduce the frequency 
difference Qj = ujj — lus- Using Af2 = 2nc/L, the transi- 
tion to the continuum limit is performed via the replace- 
ment J2j (Af2) -1 / 2 dfia(w5 + Q) ... where 
the continuous field operators a(w) have the dimension 
s 1 / 2 . Since the bandwidth of the signal is small compared 
to its central frequency lus, the integration interval can 
be extended to — oo, and we arrive at the approximate 
operator representation 



4 +3 om) = 



where 125 



2e cA 



V5n 



[a(wi),a t (w2)] = S(ui - io 2 ). 



(2.5) 



(2.6) 



Similarly, the corresponding operator for the idler field 
in free space is given by 



E { j + \x,t) 



where 



hioi 
2e Q cA 



b(uji + n)c l(uJi+n ^- t \ 



[b(LJi),b\u) 2 )} = S{lji - u 2 ). 



(2.7) 



(2. 



Since in type-II parametric down-conversion signal and 
idler photons are polarized orthogonally, we have 

[o(w 1 ),6t( Wa )]=o. (2.9) 



C. The field operators inside the crystal 

We now turn to the fields inside the crystal. If the non- 
linear interaction is small, we do not need to consider the 
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complicated problem of field quantization in a nonlinear 
medium, but we can represent the field operators inside 
the crystal by adapting the corresponding expressions for 
the free-field operators in order to account for the pres- 
ence of a lossless dispersive medium In accordance 
with Ref. 24], the operator for the positive- frequency 
part of the idler field inside the crystal then takes the 
approximate form 



2eocAni J _ 



dfl 
V5n 



(2.10) 

j-^ e i[fci(0)a;-(w.f+f2)i] 



where we introduced the wave vector at frequency luj + 91, 



9. 



(2.11) 



Here, the replacements eo ~~ * e o r *i and c — * c/n^ have 
been performed where rii is the refractive index of the 
idler wave and ni = ni{u>i). 

To describe the signal field, we have to take the pres- 
ence of the resonator into account. Let us first assume 
a lossless resonator completely filled with the nonlinear 
medium. The quantization length of the field is then 
equal to the crystal length I. When n s denotes the re- 
fractive index of the signal, the adapted resonator eigen- 
frequencies characterizing the longitudinal modes can be 
written as 

(m + m)nc UJ S n s l , . 

w m = n— with m = , (2.12) 

n s {uj m ) I ire 

where ns = n s (u>s), m = 0, ±1,±2, ... and mo ^> |m|. 
Since the frequency difference between adjacent modes is 
small compared to the total spectral width of the signal, 
as will become obvious in Sec. IV, we can assume without 
lack of generality that wg coincides with the frequency of 
a longitudinal mode and mo is an integer, i.e. that the 
cavity is tuned to resonance. Using the Taylor expan- 

sion nsjoJ m) = n s + (u> m - ^s)^f \u=u s > we nnd from 
Eq. (|2.12p after minor algebra that 



us 



7r%s 



m Au) r 



(2.13) 



where 



IJ g,S 



ns 



is the group velocity of the signal at frequency u>s- By 
adapting the empty-cavity field operator [23[ with the 



replacements eo 



and 



c/n s , the part of 



the standing-wave signal field operator inside the crys- 
tal that corresponds to a component traveling in positive 
x-direction in the lossless resonator is found to be 



4 + ci(M) 



hujsAuj c 
eonscAir 



3 iw m [f n a {u m )-t)] 



m — — oo 



(2.15) 



Here, we replaced the quantization length under the 
square-root sign by the expression I — tt v 9i s / Aw c ~ 
TTc/(n s Aoj c ), following Eqs. (|2~T3]) and (|2~14l) . Moreover, 
the summation has been extended to m = — oo, in anal- 
ogy to the expanded integration range in Eq. (|2.5p . The 
photon annihilation and creation operators for mode m 
obey the usual commutation relation [a m ,o' m i] — 5 m , m '- 
When resonator losses are incorporated, the modes 
turn into quasi-modes and the annihilation operators a m 
in Eq. (|2.15p become time-dependent . According to the 
input-output formalism [23l . |25| for a one-sided cavity 
with loss constant 7, the damping of mode m is described 
by 



a m {t) 



7 



a m (t) + ^a^{t) 



(2.16) 



The operators a 1 ^ (t) and n„ (i) characterize the ingo- 
ing and outgoing photon flux at frequency uj m and have 
the dimension s -1 / 2 . They are related by the boundary 
condition 

a™(t)=V7am(t)-ag UT (t). (2.17) 
In order to determine a m (t), we use the representation 
1 



,OUT 



(*) = 



dQ a(oj m + il)c~ 



iQI 



(2.18) 



where, in analogy to Eq. 

[a{aj m + n),a\u m , + n')] = *m,m'<y(n - n'). (2.19) 



Eq. (|2.19p implies that the quasi-modes do not overlap 
which is justified in the good-cavity limit 



7< Aw c 



(2.20) 



After inserting Eq. (|2.17p into Eq. (|2.16p . we obtain by 
Fourier transformation the solution 



a m (t) 



1 



7 27T J- 



dfl a(ujn 



'2 + in 



(2.21) 



which has to be applied to Eq. (|2.15p . The operator for 
the relevant field component of the signal in the lossy 
cavity can then be written as 



(2.14) 4+ } r (x,t) 



ftw s ^lAu c 



E 



2eonscA 2tt 

a{io m + fl) i[k s m (n)x-{w m +n)i\ 

1 + 



(2.22) 



dQ 



Here. 



(fi) 



n s (w m + Q) 



(2.23) 



is the wave vector corresponding to a traveling-wave com- 
ponent of frequency ui m + 91. The denominator in the 
integral in Eq. (|2.22p describes radiation suppression for 
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frequencies u> with \u> — u) m \ = |f2| 3> 7 while resonance 
enhancement occurs for lu ss uj m . 

So far, we have assumed a resonator length L r that 
coincides with the crystal length I. If L r > I, a rigor- 
ous quantization of the signal field has to account for the 
exact position of the crystal inside the resonator, but is 
beyond the scope of the present paper. For the purposes 
of our approximative treatment, however, it is sufficient 
to describe the signal field inside the crystal by Eq. (|2.22|) 
with Auj c = Aw and ui m — Lus+mAu; where Aw is the ef- 
fective free spectral range. The latter can be represented 



2tt . , m 21 2(L r - 1) 
Alu = — with T = h -, (2.24) 

T Vg jS c 

where T is the effective cavity round-trip time of a signal 
photon. 



III. THE RATE OF BIPHOTON GENERATION 
AND THE BIPHOTON WAVE FUNCTION 

With the expressions for the operators of the signal and 
idler field at hand, we are now in the position to specify 
the interaction Hamiltonian and to derive a perturba- 
tive solution of the Schrodingcr equation. Making use of 
Eqs. (EU), (TO]), (T2~lTJ|) and (|2~2^|) . as well as Eq. (|2~T51) 
with Atu c — > Alo and taking the frequency-dependence 
of the nonlinear susceptibility into account, we find from 

Eq. mm 



Hi, 



m— — 00 



dn' F m (n,n') 



x a){uo m + fl)tf(uji + n')e l( - mAuJ+n+n,)t + H.A., (3.1) 
where we introduced the function 

(3.2) 



F m (n,n 7 ) = 



x(w P ;w s ,w ; ) I 
and defined the constant 

-iEp ju>su>i 



d x e i[iip-V(")-*iW 



8neocA \j nsnj 



x(wp; uj s , wj)vAw. 



(3.3) 



We are interested in the regime far below threshold where 
the biphoton production rate k is much smaller than the 
cavity damping rate, 



« < 1, 



(3.4) 



and the mean photon number in the resonator there- 
fore close to zero. Since the mean time interval between 
biphoton emission events is large compared to the cav- 
ity damping time, the resonator can be assumed to be 
empty before each emission event. For this case, we can 



perform a perturbative treatment of the nonlinear inter- 
action, assuming that at the initial time t = 0, the com- 
bined signal-idler field is in the vacuum state, described 
by |0) = \0)s 8> |0)j. In the following we rely on the 
ideas developed for the th eory of spontaneous paramet- 
ric down-conversion [l], l22j |. By expanding the formal 
solution of the time-dependent Schrodinger equation, up 
to a normalization factor, the state vector describing the 
combined signal- idler field at a time t = At <C k^ 1 is 
found to be \ip(At)} cx |0) + |^(At)) with 



1 



At 



|^(At)) = - / dtH mt (t)\Q). 



(3.5) 



After substituting Eq. ()3.1j) into Eq. (|3.5p . integration 
with respect to t yields 

00 />oo rr- />oc 



i(mAw + tt + n')At 



x At sine 
xa t (w m + O)6 t (a; z + n')|0), 



e j(mAu+0+!l')At 

(3.6) 



where we used the sinc-function defined as sinc(z) = 
with sinc(0) = 1. 

The nonnormalized vector \ijj(Ai)) refers to the state 
of the radiation field on the condition that a signal-idler 
photon pair has been produced during the time inter- 
val At. Since this condition applies with the probability 
(ip(At)\ip(At)) , the biphoton production rate is given by 



K= s (V(A<MAi)). 

Applying the commutation relations Eqs. 
(pn9|) . we get from Eq. (HU) 



(3.7) 
and 



00 ^>oo 

$(A*)$(Ai)) = (At) 2 / dn 



l«| 2 7 



m— — 00 



(i) 2 + fi 2 



(3.8) 



/CO 
dCl' \F m (n,n')\ 2 sine 5 
-OO 



-(mAuj + tt + n')At 



Because of the properties of the sinc-function, the inte- 
gral is dominated by the region mAuj + il + fl' < -k j At. If 
At is sufficiently large and F m (fl, f2') is a slowly varying 
function of O' in this region, the latter function can be 
replaced by its value at ft' — -mAw — Q 1]. Using the 
relation sinc 2 (z)iiz = tt, the integration with respect 
to Q' is then readily performed, and we obtain 



$(At)WAt)) = 
2n\a\ 2 At ]T / 



(3.9) 



dfl 



(i) 2 



n 2 



where 



$ m (fi) = F m (Cl, -mAw - Q). 



(3.10) 
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For further evaluation, we need to specify the func- 
tion $ ro (0). Using Eqs. $FZ\), (gXp , and (|2~2l?)) with 
u.,„ = lus + mAw, Taylor expansion around the central 
frequencies u>s and u>i yields 

kp - fcs, ro (0) - fc/(-mAw-0) w (mAw + 0)y, (3.11) 
where we introduced the time constant 



m + uji 



drij 



n s - cos 



dn s 
du> 



The latter is equivalent to 



I 



TO = 



V 9,I V 9,S 



(3.12) 



(3.13) 



and describes the difference between the transit times 
of a signal and idler photon through a crystal of length 
I, originating from the difference in the signal and idler 
group velocities, v g g and v g j, defined by Eq. (|2.14p and 
by the corresponding equation for the idler wave, respec- 
tively. Since in any real experiment \v 9i s — v g j\ <C v 9l s, 
it follows that |t | < (Aw) -1 where we used Eq. (|2.24|) . 
Hence, we are considering a parameter range in this prob- 
lem that is characterized by the combined inequality 



k < 7 < Aw« |t | 



(3.14) 



where Eqs. (|2.20p and (|3.4p have been incorporated [27j |. 
Assuming a constant nonlinear susceptibility within the 
bandwidth given by |r | _1 , we find from Eqs. (|3.10[) . 
([3~TTj) . and that 



1 f° 
$m(ty™jj dx 



i(mAu+!))3.i 



(3.15) 



In order to determine the rate of biphoton generation, 
we have to insert Eq. (|3.15p into Eq. (|3.9[) where the in- 
tegration with respect to is effectively restricted to the 
cavity bandwidth 7 with 7 <C Aw. Hence, we neglect 
f2 compared to mAw in Eq. (|3.15[) for m 7^ 0. More- 
over, Eq. (|3.14p implies \£It \ <C 1 in the relevant interval 
|Q| < 7 and therefore $o(^) ~ 1- After integration with 
respect to we get the approximation 



$ m (ft) a $ m (0) = sine ( mA U |) e ~*™Aui 



(3.16) 



which can be used in connection with Eq. (|3.9p , (|3.7p . 
and (|3 . 3|) to determine the rate of biphoton generation 



Aw y. smc 



( mAw— ] 

V 2 / 



4encA / riou/ 

v / m— — 00 

(3.17) 

Further simplification is possible if we transform the sum 
into an integral with respect to z = mAwro/2, introduc- 
ing the positive increment dz — Ac<j|to|/2 where dz <C 1 
because of Eq. (|3.14| . We then arrive at the expression 



\4e c^4 / nsni \tq\ 



(3.18) 



that does not depend on the properties of the resonator 
because the mean photon number in our cavity is approx- 
imately zero and, in addition, an associated idler mode 
exists for each signal mode which meets the requirement 
for energy conservation. 

The presented perturbative treatment relies on the 
condition |to| <C At <C k _1 so that the wave function 
l^i(At)) in Eq. (|3.6p and the approximation leading to 
Eq. (|3.9[) are valid simultaneously. The normalized vector 
= {K,At)^ 1 / 2 \ijj{At)) can be denoted as the biphoton 
wave function since it represents the state of the radia- 
tion field on the condition that exactly one signal-idler 
photon pair is present. 

A more convenient representation of the biphoton wave 
function is obtained if the sinc-function in Eq. (|3.6p is 
replaced by 2ir5(mAuj+Q+Sl.') according to the standard 
procedure Mathematically, this corresponds to the 
limit At — > 00, implying k — > 0. Then, the integration 
with respect to fi' can be performed immediately. In 
analogy to the expression given in Ref. 22|, we obtain 
the biphoton wave function 



E / 

„ j — 1 



dn 



gm(fl) 



+ n)tf{u I -mAn-n)\o), (3.19) 



where $ m (S!) is given by Eq. (|3 . 1 5[) and Af is a normaliza- 
tion constant. The explicit value of Af (2(| is not impor- 
tant as long as only normalized quantities characterizing 
the radiation field are considered. Eq. (|3.19p clearly re- 
veals the frequency-entanglement between the signal and 
idler photon and will serve as our basic equation to de- 
termine the properties of the emitted radiation. 



IV. PROPERTIES OF THE EMITTED 
RADIATION 

A. Output spectra of the signal and idler field 

First we investigate the output spectra Ss(to) and 
Si(uj) of the signal and idler field, defined as 



(i) 

s/i 



(r) e* 



(4.1) 



where G^^r) and G^p (r) are the first-order temporal 
correlation functions of the respective fields outside the 
resonator. In the following, it is convenient to use the 
Heisenberg picture and to start from the expression 

G«(r) = (V>|47](x,t)4+>(z,i + T)|V>>. (4.2) 
Here, \ip) is the time-independent biphoton wave func- 



tion given by Eq. (|3~Ti?j) and E ( s +> = E ( s 11 and E { j +) = 

(-)+ 

Ej denote the positive-frequency parts of the time- 
dependent operators of the electric fields in free space, 
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-2*/|i | 

Frequency co-co s/ | 



2«/|t l 



at frequencies u>s/i ~ mAto with m = 0, ±1, . . . where 
the respective spectral envelopes are determined by the 
sinc-function in the nominator. From Fig. 2 it becomes 
obviouos that the frequency bandwidth of the signal and 
idler photons is characterized by |to| _1 - Their tempo- 
ral uncertainty is thus equal to the modulus of the time 
constant ro introduced in Eq. (|3.12|) and resulting from 
the phase-matching conditions. Even though the idler 
wave is not resonant, the longitudinal mode structure of 
the resonator is revealed in the idler spectrum due to 
the frequency entanglement between the signal and idler 
photon which arises from the interaction underlying the 
biphoton generation process. 



B. Signal-idler cross-correlations 



FIG. 2: (Color online) Schematic plot of the normalized sig- 
nal and idler output spectra. The width of the Lorentzian 
peaks is determined by the cavity damping rate 7, and they 
are separated by the free spectral range Auu = 2n/T. The 
envelope, described by the sinc-function in Eq. (|4.6f) . yields 
the total spectral width 2tt/\tq\. 



given by Eqs. (|2.5p and (|2.7[) . Let us first determine the 
spectrum of the idler field. Considering 

b(uJi + ft')^/ - mAw - O)|0) = S(fl' + mAuo + Q)\0), 

(4.3) 

due to Eq. 1^8]) . we obtain from Eqs. ([2J]) and (|3~T9)) 



4 +) 0MM« £ J 



m— — 00 



Sm(fl) 



(4.4) 



x e l( ^- mA "- 0) (t-*) a t( Ws + m Auj + (l)\0). 

By taking the inner product of Ej(x,t + T)\ip) and 

Ej (x, t)\ip) and applying the commutation relation 
Eq. (|2~P9"j) . we find 



J — ( 



— — mAw- Q)t 



Q 2 



(4.5) 



In analogy to the derivation of Eq. (|3.17|) , the dependence 
of $ m on can be neglected within the relevant band- 
width determined by 7. After inserting Eq. p,16[) into 
Eq. (|4.5p . the Fourier transform yielding the idler spec- 
trum according to Eq. (|4.1|) is readily performed. The 
signal spectrum can be determined in a completely anal- 
ogous way, and we finally arrive at the relation 



Ss/i(u) oc 



sine 2 (mAw^ 1 ) 



(J) + {^s/i - mAuj - uj) 2 



(4.6) 



Eq. (|4.6|) indicates that both the signal and idler spec- 
trum are composed of Lorentzians of halfwidth 7 centered 



The coincidence rate for detecting an idler photon at 
time t and a signal photon at time t + r, both at equal 
distance from the end facet of the crystal, is proportional 
to the temporal correlation function 

Gfj(r)= (4.7) 

{^\E { f ] (x, (M + T )4 +) t + T ) E< i +) MM, 

where again E^ and Ej are the free-field operators 



defined by Eqs. (f2~5j) and (|2~7| . With the explicit ex- 
pression for the biphoton wave function given by 
Eq. dnHU), we get 



E^\x,t + T)E^ +) (x,t)\ip) cx e 4^+^)(t-t)-r WS ] 

x g c -™a- d nm^e- inr |0), (4.M 



rn— — 00 



where Eq. (|4.3[) and the corresponding relation for the 
signal modes 

a(w s + f7')a t (cj s + TOAw + ^)|0) = S(fl' -mAuj-n)\0), 

(4-9) 

following from Eq. (|2.6p . have been applied. According 



— — (2) 
to Eq. (|4.7p . the correlation function Gj S (r) is propor- 
tional to the squared norm of the Hilbert vector on the 
right-hand side of Eq. (|4.8p . It can be determined from 
the explicit expression for $ m (0), given by Eq. (|3.15p . 
together with the integral identities 



if t < 

if t = (4.10) 
if t > 0. 




First, it is important to observe that 



G$(r)=0 



if 



Hi 

2 



< 



M 
2 



(4.11) 



which is equivalent to r < —To for tq > and r < for 
r < 0, respectively. Mathematically, Eq. (|4.1ip is due to 
the expression fit = 0(r— fio) which results from insert- 
ing Eq. (|3.15p into (|4.8[) and is negative in the given case 
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Time Delay - 



FIG. 3: Schematic representation of the normalized second- 
order signal-idler cross-correlation function (r) according 
to Eq. (|4~13j) for r > and -y/Aco = 0.05. If the time delay 
t is equal to zero or multiples of the cavity round-trip time 
T = 27t/Aoj, the function exhibits pronounced peaks which 
decay with the cavity damping time 7 _1 . 



of the signal wave packet to travel from the middle of the 
crystal to its end facet is by the amount to/2 shorter than 
the time needed by the center of the idler wave packet 
to cover the same distance while the opposite holds for 
r < 0. 

Transforming the sum in Eq. (14. 1 2[) into an inte- 
gral, the expression for the signal-idler cross-correlation 
function can be further approximated, in analogy to 
the procedure applied to derive Eq. (|3.18[) . Introduc- 
ing z = mAuJTo/2 and dz = Aujtq/2, we find that 
Gjg (t) oc e _7T / 2 with I = sinc(z) cos(az) dz since 
sinc(z) is an even function of z. Because of the period- 
icity of the cos-function, the parameter a can be written 
as a = 1 + 2 [r -([£]+ 1) T] /r where [ denotes the 
largest integer that does not exceed t/T. Since I = tt for 
\a\ < 1 and I = for \a\ > 1, Eq. (|4.12j) takes a simple 
form that can be combined with Eq. (|4. 1 1[) to yield the 
compact approximate representation 



G§(t)c< 



oo 

E 

3=0 



-7jT 



if | 

else, 



Id. 

2 



< 



M 

2 



(4.13) 



for any x inside the crystal, i.e. for — I < x < 0. There- 
fore, the upper line of Eq. (|4.10p applies. Physically, the 
correlation function vanishes because the arrival time of 
the signal photon in a photon pair can precede the arrival 
time of the corresponding idler photon at most by a time 
interval that is within the temporal uncertainty interval 
| To | inherent in the biphoton generation process. 

On the other hand, the presence of the resonator al- 
lows a signal photon to be detected considerably later 
than the associated idler photon since the signal photon 
may bounce back and forth between the resonator mir- 
rors repeatedly before leaving the resonator. For time 
delays r outside the limits of Eq. (|4.11[) . we approxi- 
mate & m (fl) by Eq. (|3.16[) where the integration with 
respect to x has been performed and the dependence 
on £1 has been neglected. After inserting Eq. (|3.16p 
into Eq. (|4.8[) . the damping of the cavity gives rise to 
a factor proportional to e~ ? T for Gis(t) j^z Q which 
follows from the third line of Eq. (f4~T0)) (28j . Since 



from Eq. 



f\x,t + T)E^ + \x,t)m 

the approximate result 



Gf s {r) ex e"^ 



if 



oc 

sine [riiAui-^- 



we finally get 



m— — oo 



Jo. > _ 
2 — 



Jjol 
2 ' 



(4.12) 



According to Eq. (]3.13[) . the sign of r is determined by 
the relation between the signal and idler group velocities, 
i.e. t > for v 9: s > v g j and t < for v g ,s < v g j. A 
numerical evaluation reveals that Eq. (|4.12[) describes a 
decaying periodic function with peaks of width |to | cen- 
tered at t — j T—t /2 where j = 0, 1, . . . and T = 2tt/ Auj 
is the cavity round-trip time. The time shift —tq/2 of the 
peaks arises since for tq > the time needed by the center 



where j = | [jr] | (see Fig. 3). Here, we took into account 
that due to Eq. (|3.14p cavity damping is negligible during 
the time interval |tq | , i.e. exp(— 7T0) ~ 1. 



DISCUSSION AND CONCLUSIONS 



In a real experiment, the sharp peaks of the function 

( 2) 

Gjg(r) will be broadened due to the finite resolution 
time of the detector setup. When the latter is taken into 
account by performing the convolution with respect to 
a Gaussian function, Eq. (14.13|) yields the time-averaged 




2 4 

Time Delay t/T 



FIG. 4: Time-averaged second-order signal-idler cross-corre- 

(2) 

lation function Gis (t) for 7/Acj = 0.05. The ratio be- 
tween the resolution time AT and the cavity round-trip time 
T = 2n/Aijj is assumed as AT /T = 0.02 (dashed line) and 
AT /T = 1 (solid line), respectively. 
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cross-correlation function 



G 



IS 



-(2) 



(r) oc 



E 

j=0 



exp 



4(jT- 



^21 



(AT)- 



(5.1) 



where AT characterizes the effective resolution time and 
where we have assumed AT 3> |to|. The resulting aver- 
aged function is plotted in Fig. 4 for two different values 
of AT. A second-order cross-correlation function show- 
ing the behavior of the solid line in Fig. 4 has recently 
been measured in our group (20| . and the results have 
been found to be in excellent agreement with the predic- 
tions derived from Eqs. (|4.1ip and (|4.12jl . 

We still note that for spontaneous parametric down- 
conversion in a double-resonant cavity the signal-idler 
cross-correlation function has also been found to exhibit 
a comb-like structure which is. however, symmetric with 
respect to the time delay Q EH, [IB]. The effect has 
been explained by applying the concept of mode-locking 
to the frequency-entangled biphoton state, pointing out 
that due to the large coherence time of the pump, pho- 



ton pairs with different frequencies have a common phase 
and form a coherent superposition [Til ]. 

To summarize, we performed a theoretical investiga- 
tion of biphoton generation by spontaneous paramet- 
ric down-conversion in a single-resonant OPO far below 
threshold. We derived analytical expressions for the rate 
of biphoton generation, for the output spectra of the sig- 
nal and idler fields, as well as for the second order signal- 
idler cross correlation function. Our investigations pro- 
vide the theoretical background for explaining the results 
of a recent experiment [201 ] , where stable continuous op- 
eration of a single-resonant OPO far below threshold has 
been demonstrated. 
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